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INTRODUCTION 



Topological objects, in particular finite energy topo- 
logical objects, have played important roles in physics 
Pi 0] ■ In Bose-Einstein condensates (BEC) the best 
known topological objects are the vortices, which have 
been widely studied in the literature. Theoretically these 
vortices have successfully been described by the Gross- 
Pitaevskii Lagrangian. On the other hand, the recent ad- 
vent of multi-component BEC, in particular the spin-1/2 
condensate of 87 Rb atoms, has widely opened an interest- 
ing possibility for us to construct totally new topological 
objects in condensed matter physics. This is because 
the multi-component BEC obviously has more interest- 
ing non-Abelian structure which does not exist in ordi- 
nary (one-component) BEC, and thus could admit new 
topological objects which are absent in ordinary BEC 
0, 0] . As importantly, the multi-component BEC pro- 
vides a rare opportunity to study the dynamics of the 
topological objects theoretically. The dynamics of multi- 
component BEC could be significantly different from that 
of ordinary BEC. This is because the velocity field of the 
multi-component BEC, unlike the ordinary BEC, in gen- 
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eral has a non-vanishing vorticity which could play an 
important role in the dynamics of the multi-component 
BEC So the multi-component BEC provides an ex- 
cellent opportunity for us to study non-Abelian dynamics 
of the condensate theoretically and experimentally. 

The purpose of this paper is to discuss the non- 
Abelian dynamics of two-component BEC. We first study 
the popular Gross-Pitaevskii theory of two-component 
BEC, and compare the theory with the recent gauge the- 
ory of two-component BEC which has a vorticity interac- 
tion §. We show that, in spite of the obvious dynamical 
differences, two theories are not much different physi- 
cally. In particular, they admit remarkably similar topo- 
logical objects, the helical vortex whose topology is fixed 
by Ti2{S 2 ) and the vorticity knot whose topology is fixed 
by tt^(S 2 ). Moreover we show that the vorticity knot is 
nothing but the vortex ring made of the helical vortex. 
Finally we show that the gauge theory of two-component 
BEC is very similar to the theory of two-gap supercon- 
ductors, which implies that our analysis here can have an 
important implication in two-gap superconductors. 

A prototype non-Abelian knot is the Faddeev-Niemi 
knot in Skyrme theory |(| @] . The vorticity knot in two- 
component BEC turns out to be surprisingly similar to 
the Faddeev-Niemi knot. So it is important for us to 
understand the Faddeev-Niemi knot first. The Faddeev- 
Niemi knot is described by a non-linear sigma field h 
(with h 2 = 1) which defines the Hopf mapping ^(S 2 ), 



2 



the mapping from the compactified space S 3 to the target 
space S 2 of n, in which the preimage of any point in the 
target space becomes a closed ring in S 3 . When ^(S 2 ) 
becomes non-trivial, the preimages of any two points in 
the target space are linked, with the linking number fixed 
by the third homotopy of the Hopf mapping. In this case 
the mapping is said to describe a knot, with the knot 
quantum number identified by the linking number of two 
rings. And it is this Hopf mapping that describes the 
topology of the Faddeev-Niemi knot 



In this paper we show that the vorticity knot in two- 
component BEC has exactly the same topology as the 
Faddeev-Niemi knot. The only difference is that here 
the vorticity knot in two-component BEC has the extra 
dressing of the scalar field which represents the density 
of the condensation. 

The paper is organized as follows. In Section II we 
review the Skyrme theory to emphasize its relevance in 
condensed matter physics. In Section III we review the 
topological objects in Skyrme theory in order to com- 
pare them with those in two-component BEC. In Sec- 
tion IV we review the popular Gross-Pitaevskii theory 
of two-component BEC, and show that the theory is 
closely related to Skyrme theory. In Section V we dis- 
cuss the helical vortex in Gross-Pitaevskii theory of two- 
component BEC, and show that it is a twisted vorticity 
flux. In Section VI we discuss the gauge theory of two- 
component BEC which includes the vorticity interaction, 
and compare it with the Gross-Pitaevskii theory of two- 
component BEC. In Section VII we discuss the helical 
vortex in gauge theory of two-component BEC, and com- 
pare it with those in Gross-Pitaevskii theory and Skyrme 
theory. We demonstrate that the helical vortex in all 
three theories are remarkably similar to one another. In 
Section VIII we present a numerical knot solution in the 
gauge theory of two-component BEC, and show that it is 
nothing but the vortex ring made of helical vorticity flux. 
Finally in Section IX we discuss the physical implications 
of our result. In particular we emphasize the similarity 
between the gauge theory of two-component BEC and 
the theory of two-gap superconductor. 



II. SKYRME THEORY: A REVIEW 



The Skyrme theory has long been interpreted as an 
effective field theory of strong interaction with a re- 
markable success |9j- However, it can also be inter- 
preted as a theory of monopoles, in which the monopole- 
antimonopole pairs are confined through a built-in Meiss- 
ner effect 0, |g . This suggests that the Skyrme theory 
could be viewed to describe a very interesting condensed 
matter physics. Indeed the theory and the theory of two- 
component BEC have many common features. In partic- 
ular, the topological objects that we discuss here are very 
similar to those in Skyrme theory. To understand this we 



review the Skyrme theory first. 

Let lo and n (with n 2 — 1) be the massless scalar field 
and non-linear sigma field in Skyrme theory, and let 



U = cxp(— a ■ n) — cos — — i(a • n) sin — , 
2i 2 2 
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With this one can write the Skyrme Lagrangian as Q 
C 



-E 
4 

a 
16 



^tr Ll + gtv([L^L v ]f 

i(M 2 + 2sin 2 |(9,n) 2 
sin 2 —(dpLodvh — d v Lod^n) 2 

2^ 7 



+4 sin 4 — (dfifi x d u n) 



(2) 



where fj, and a are the coupling constants. The La- 
grangian has a hidden local E/(l) symmetry as well as 
a global SU (2) symmetry. From the Lagrangian one has 
the following equations of motion 
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Notice that the second equation can be interpreted as a 
conservation of SU (2) current, which of course is a simple 
consequence of the global SU (2) symmetry of the theory. 

With the spherically symmetric ansatz 

lo = Lo(r), h = f, (4) 

(3) is reduced to 
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Imposing the boundary condition 

lo(0) = 2tt, w(oo) = 0, 



(5) 
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one can solve the Eq. (JSJ and obtain the well-known 
skyrmion which has a finite energy. The energy of the 
skyrmion is given by 
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where a; = u/^fa is a dimensionless variable. Further- 
more, it carries the baryon number (3-0 
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which represents the non-trivial homotopy 713 (iS* 3 ) of the 
mapping from the compactified space S 3 to the SU(2) 
space S 3 defined by U in Q. 

A remarkable point of © is that 



Cl) = 71", 



(9) 



becomes a classical solution, independent of n [(J. So 
restricting cj to 7r, one can reduce the Skyrme Lagrangian 
to the Skyrme-Faddeev Lagrangian 
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whose equation of motion is given by 
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Notice that -ff^j, admits a potential C M because it forms 
a closed two-form. Again the equation can be viewed as 
a conservation of SU(2) current, 



df_,(n x 9 M n + —H^ v d v n) = 0. 
A* 



(12) 



It is this equation that allows not only the baby skyrmion 
and the Faddeev-Niemi knot but also the non-Abelian 
monopole |6j,|8(. 



III. TOPOLOGICAL OBJECTS IN SKYRME 
THEORY 



The Lagrangian (|10(l has non-Abelian monopole solu- 
tions lal 



n = r, 



(13) 



where r is the unit radial vector. This becomes a solution 
of l|ll|) except at the origin, because 



d 2 r = -^r, dnHuu = 0. 



(14) 



20 



41) 



60 



80 



100 



FIG. 1: The baby skyrmion (dashed line) with m — 0,n — 1 
and the helical baby skyrmion (solid line) with m = n = 1 in 
Skyrme theory. Here q is in the unit y/a/ ^ and k = 0.8 \xj \fa. 



This is very similar to the well-known Wu-Yang 
monopole in SU(2) QCD 0, It has the magnetic 

charge 
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which represents the non-trivial homotopy ^(S 2 ) of the 
mapping from the unit sphere S 2 centered at the origin 
in space to the target space S 2 . 

The above exercise tells that we can identify H^ v a 
magnetic field and C M the corresponding magnetic po- 
tential. As importantly this tells that the skyrmion is 
nothing but a monopole dressed by the scalar field u>, 
which makes the energy of the skyrmion finite |(J . 

It has been well-known that the Skyrme theory has a 
vortex solution known as the baby skyrmion \\\\ . More- 
over, the theory also has a twisted vortex solution, the 
helical baby skyrmion Q . To construct the desired heli- 
cal vortex let (g, ip, z) be the cylindrical coordinates, and 
choose the ansatz 



(sin f(g) cos (nip + mkz) 1 
sin f(g) sin (nip + mkz) 
cos /(g) 



(16) 



With this we have (up to a gauge transformation) 

C M = -(cos/ + l)(nd^ip + mkdfj.z), (17) 
and can reduce the equation 111|) to 
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FIG. 2: The supercurrent i$ (in one period section in z- 
coordinate) and corresponding magnetic field H z circulating 
around the cylinder of radius g of the helical baby skyrmion 
with m — n — 1. Here g is in the unit \/a/fi and k = 
0.8 m 'yfa. The current density j$ is represented by the dot- 
ted line. 



FIG. 3: The supercurrent i$ and corresponding magnetic field 
H<p flowing through the disk of radius g of the helical baby 
skyrmion with m = n = 1. Here g is in the unit \fa/iJt and 
k — 0.8 (x/^/a.. The current density j z is represented by the 
dotted line. 



So with the boundary condition 

/(0)=7T, /M=o, 



(19) 



we obtain the non-Abelian vortex solutions shown in 
Fig. ^ Notice that, when m = 0, the solution de- 
scribes the well-known baby skyrmion. But when m is 
not zero, it describes a helical vortex which is periodic 
in z-coordinate Q- In this case, the vortex has a non- 
vanishing magnetic potential C M not only around the vor- 
tex but also along the z-axis. 

Obviously the helical vortex has the helical magnetic 
field made of 

1 n ■ 

H z = -H ev = -/sin/, 
g g 
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which gives two quantized magnetic fluxes. It has a quan- 
tized magnetic flux along the z-axis 



= / H ev dgdip — —Ann, 



(21) 



and a quantized magnetic flux around the z-axis (in one 
period section from to 27r/fc in z-coordinate) 



= — I H gz dgdz = 4irm. 



(22) 



Furthermore they are linked since </>i is surrounded by 
4><p. This point will be very important later when we 
discuss the knot. 

The vortex solutions implies the existence of Meissner 
effect in Skyrme theory which confines the magnetic flux 
of the vortex To see how the Meissner effect comes 



about, notice that due to the U(l) gauge symmetry the 
Skyrme theory has a conserved current, 



j M = dyH^, d^jf, = 0. 



(23) 



So the magnetic flux of the vortex can be thought to 
come from the helical electric current density 



^ = -sin/L(/ + ^/ 2 -i/)^ 
L v sin/ Q 

sin/ Q 



(24) 



This produces the currents (per one period section in 
z-coordinate from z = to z = 2n/k) around the z-axis 
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and iz along the z-axis 
it = — mk 



e=co ■■ cos/ • 1 ■ 

sin/(/ + -^—fJ 2 + -f)gdgdif 
sin/ g ' 



lg=0 

= —2irmkgfsmf 



= 0. 



(26) 



£>=0 



Notice that, even though iz — 0, it has a non-trivial 
current density which generates the net flux tp^. 

The helical magnetic fields and currents are shown in 
Fig. |21 and Fig- EI Clearly the helical magnetic fields are 
confined along the z-axis, confined by the helical current. 
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This is nothing but the Meissner effect, which confirms 
that the Skyrme theory has a built-in mechanism for the 
Meissner effect. 

The helical vortex will become unstable and decay to 
the untwisted baby skyrmion unless the periodicity con- 
dition is enforced by hand. In this sense it can be viewed 
as unphysical. But for our purpose it plays a very im- 
portant role, because it guarantees the existence of the 
Faddeev-Niemi knot in Skyrme theory 0, 0] ■ This is be- 
cause we can naturally enforce the periodicity condition 
of the helical vortex making it a vortex ring by smoothly 
bending and connecting two periodic ends together. In 
this case the periodicity condition is automatically im- 
plemented, and the vortex ring becomes a stable knot. 

The knot topology is described by the non-linear 
sigma field h, which defines the Hopf mapping from the 
compactified space S 3 to the target space S 2 . When the 
preimages of two points of the target space are linked, 
the mapping tt^ (S 2 ) becomes non-trivial. In this case the 
knot quantum number of ^(S* 2 ) is given by the Chern- 
Simon index of the magnetic potential C M , 



Qu 



l 



32tt 2 



(27) 



Notice that the knot quantum number can also be under- 
stood as the linking number of two magnetic fluxes of the 
vortex ring. This is because the vortex ring carries two 
magnetic fluxes linked together, m unit of flux passing 
through the disk of the ring and n unit of flux passing 
along the ring, whose linking number becomes ran. This 
linking number is described by the Chern-Simon index of 
the magnetic potential Q- 

The knot has both topological and dynamical stabil- 
ity. Obviously the knot has a topological stability, be- 
cause two flux rings linked together can not be discon- 
nected by any smooth deformation of the field. 

The dynamical stability follows from the fact that the 
supercurrent 124(1 has two components, the one moving 
along the knot and the other moving around the knot 
tube. Clearly the current moving along the knot gener- 
ates an angular momentum around the z-axis which pro- 
vides the centrifugal force preventing the vortex ring to 
collapse. Put it differently, the current generates the m 
unit of the magnetic flux trapped in the knot disk which 
can not be squeezed out. And clearly, this flux provides 
a stabilizing repulsive force which prevent the collapse of 
the knot. This is how the knot acquires the dynamical 
stability. It is this remarkable interplay between topology 
and dynamics which assures the existence of the stable 
knot in Skyrme theory 0]. 

One could estimate the energy of the knot. Theoret- 
ically it has been shown that the knot energy has the 
following bound 



where c = 8tt 2 x 3 3 / 8 and C is an unknown constant not 
smaller than c. This suggests that the knot energy is 
proportional to Q 3 ^ 4 . Indeed numerically, one finds |l8j | 
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(29) 



up to Q — 8. What is remarkable here is the sub- linear 
Q-dependence of the energy. This means that a knot 
with large Q can not decay to knots with smaller Q. 



IV. GROSS-PITAEVSKI THEORY OF 
TWO-COMPONENT BEC: A REVIEW 



The creation of the multi-component Bose-Einstein 
condensates of atomic gases has widely opened new op- 
portunities for us to study the topological objects ex- 
perimentally which so far have been only of theoretical 
interest. This is because the multi-component BEC can 
naturally represent a non-Abelian structure, and thus can 
allow far more interesting topological objects. Already 
the vortices have successfully been created with differ- 
ent methods in two-component BECs [jj 0. But the- 
oretically the multi-component BEC has not been well- 
understood. In particular, it need to be clarified how dif- 
ferent the vortices in multi-component BEC are from the 
well-known vortices in the single-component BEC. This 
is an important issue, because the new condensates could 
have a new interaction, the vorticity interaction, which 
is absent in single-component BECs. So in the following 
we first discuss the vortex in the popular Gross-Pitaevskii 
theory of two-component BEC, and compare it with that 
in gauge theory of two-component BEC which has been 
proposed recently [j| . 

Let a complex doublet 4> — (4>i,4>2) be the two- 
component BEC, and consider the non-relativistic two- 
component Gross-Pitaevskii Lagrangian [Til fljl IbH Il7| 
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h r 



-(4(d t <h) - (dtfaVfc) - 

-Ai 2 (^i)(4<fe) - ^(4<K> 2 , ,;->„, 



2M 
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where fj,i are the chemical potentials and \j are the quar- 
tic coupling constants which are determined by the scat- 
tering lengths aij 



A; 



4irh 2 



(31) 



n Q 3/4 <e q <c 



)3/4 



(28) 



The Lagrangian (|30() is a straightforward generalization 
of the single-component Gross-Pitaevskii Lagrangian to 
the two-component BEC. Notice that here we have ne- 
glected the trapping potential. This is justified if the 
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range of the trapping potential is much larger than the 
size of topological objects we are interested in, and this 
is what we are assuming here. Clearly the Lagrangian 
has a global U(l) x U(l) symmetry. 

One could simplify the Lagrangian (|3U[1 noticing the 
fact that experimentally the scattering lengths often have 
the same value. For example, for the spin 1/2 conden- 
sate of 87 Rb atoms, all have the same value of about 
5.5 nm within 3 % or so 0, EI- In this case one may 
safely assume 



An — A12 ~ A22 — A. 
With this assumption l|3l)[) can be written as 



(32) 



r h 

C = l 2 
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H = in, 



~{d t 4> 

- H) 2 
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(33) 



(34) 



Clearly the Lagrangian has a global U(2) symmetry when 
5[i = 0. So the Sp interaction can be understood to 
be the symmetry breaking term which breaks the global 
U(2) symmetry to U(l) x U(l). Physically 5fx represents 
the difference of the chemical potentials between <\>\ and 
4>2 (Here one can always assume Sp > without loss of 
generality), so that it vanishes when the two condensates 
have the same chemical potential. Even when they differ 
the difference could be small, in which case the symme- 
try breaking interaction could be treated perturbatively. 
This tells that the theory has an approximate global U (2) 
symmetry, even in the presence of the symmetry breaking 
term |18| . This is why it allows a non- Abelian topological 
objects. 

Normalizing <j> to {y2M /%)(f) and parametrizing it by 



7T C 



7T' ctc 



i) 



(35) 



we obtain the following Hamiltonian from the Lagrangian 
(|33|l in the static limit (in the natural unit c = K = 1), 



n= 1 -(d l pf+ 1 -pMc\ 2 + ^P 2 -P 2 o) 2 

+^rP C2C2, 



(36) 



where 



A = 4M*A, pl^jf x..» 

Minimizing the Hamiltonian we have 
'A 



dp 2 = 2MSp. (37) 



fflp \diC\ 2 p = (2 (p 2 p 2 ) + Sh 2 (GC2))p, 



{(9 2 -C t 9 2 C)+ 2^(5,-^0 

+^ 2 (CIC 2 )}Ci = o, 
L p 

-v 2 (crci)}c 2 = o, 

C t 9 l (p 2 ^C)-a i (p 2 ^C t )C = 0. (38) 

The equation is closely related to the equation we 
have in Skyrme theory, although on the surface it appears 
totally different from (|1 1|) . To show this we let 

n = C+ffC, 
C„ = -2iC f ^C, (39) 

and find 

[d,nf = 4(|^C| 2 - K^CI 2 ) = 4|^C| 2 ~ C% 

h ■ (d^n x d v n) = -2i(d^d u C, ~ c^C^O 

= dpCv - dvCp = Hp,. (40) 

Notice that here is precisely the closed two-form 
which appears in |1 l|l . Moreover, from 1)39(1 we have the 
identity 



dp + —{Cfj,n -fix dfj.fi) ■ a 



C = o. 



(41) 



This identity plays an important role in non-Abelian 
gauge theory, which shows that there exists a unique 
517(2) gauge potential which parallelizes the doublet £ 
|10|. For our purpose this allows us to rewrite the equa- 
tion of the doublet ( in into a completely different 
form. 

Indeed with the above identities we can express l|38|) 
in terms of n and C M . With (|40[1 the first equation of 
(|38fl can be written as 

d 2 P -\[{d l n) 2 + C 2 ]p=(\{p 2 -p 2 ) 



+^ 2 (CIC 2 ) )p 



(42) 



Moreover, with (|41l) the second and third equations of 
(|38ll can be expressed as 



2/ 



A + B-a)( = 0, 

A = did + + i(2C 2 *C 2 - 1)5 fi 2 , 

P 

-* do 
B = n X d 2 n + 2-^—n x difi — Cidih 
P 

-{did + 2-^-C t )n + iSp 2 k, 
P 

where k = (0, 0, 1). This is equivalent to 

A + B-h = Q, 

h x B — in x (h x B) = 0, 



(43) 



(44) 
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so that (|43|l is written as 



n x d 2 n + 2^-ri x dih — Cidih 
P 

— bp 2 k x h. 



(45) 



Finally, the last equation of (|38|) is written as 

dt( P 2 C t ) = 0, (46) 

which tells that p 2 Ci is solenoidal (i.e., divergenceless) . 
So we can always replace Ci with another field Bi 

Ci = -^tijkdjBk = —-^diGij, 



Gj,j ^ijk-Bki 



(47) 



and express iUST) as 



h x d 2 h + 2^- h x dih H — ^diGijdjh 



P P 
= bp 2 k x n. 



1 



(48) 



With this (|38|l can now be written as 

d 2 P - 1 -[(d t n) 2 +C 2 ]p=(±(p 2 -pt) 



+V(C 2 *C 2 )Jp, 
h x d 2 h + 2^-h x dih + -^-diGijdjh = bp 2 ft x ft, 



P P 

% Gij 



di.Gn — —p 2 Cj. 



(49) 



This tells that (|38[) can be transformed to a completely 
different form which has a clear physical meaning. The 
last equation tells that the theory has a conserved U(l) 
current j^, 



3n = P 2 C lll 



(50) 



which is nothing but the Noether current of the global 
U(l) symmetry of the Lagrangian The second equa- 
tion tells that the theory has another partially conserved 
SU(2) Noether current j^, 



jf, = p 2 h x d^h - p 2 C fJ ,h, 



(51) 



which comes from the approximate SU (2) symmetry of 
the theory broken by the bp term. It also tells that the 
theory has one more U(l) current 



"V* — k ■ in, 



(52) 



which is conserverd even when bp is not zero. This is 
because, when Sp is not zero, the 5/7(2) symmetry is 
broken down to a U(l). 

More importantly this shows that (|38[1 is not much 
different from the equation in Skyrme theory. Indeed 
in the absence of p, Ijllfl and 149|) acquire an identical 



form when 6p 2 = 0, except that here is replaced by 
Gij . This reveals that the Gross-Pitaevskii theory of two- 
component BEC is closely related to the Skyrme theory, 
which is really remarkable. 

The Hamiltonian (|36[1 can be expressed as 



h = \pt n, 
n = l(d l p) 2 + 1 - P 2 \d l c\ 2 + l(p 2 -i) 2 



where 



+ Tp P C ^ 2 ' 



P » » 1 

P= , Oi = Kdi, K = 



(53) 



Po 



Notice that H is completely dimensionless, with only one 
dimensionless coupling constant bp I p. This tells that the 
physical unit of the Hamiltonian is \p\ , and the physical 
scale k of the coordinates is l/\f~\pQ. This is comparable 
to the correlation length £, 



£ = ; 1 = V2 K. 
^ ^/2pM 



(54) 



For 87 Rb we have 



M ~ 8.1 x 10 10 eV, A ~ 1.68 x 1CT 7 (nm) 2 , 

fj, ~ 3.3 x 10~ 12 eV, bp ~ 0.1 p, (55) 



so that the density of 87 Rb atom is given by 

< 0t0 >= ft ~ o.998 x 10 14 /cm 3 . 
A 

From (|55|1 we have 



(56) 



A~1.14xlO u , p\ ~ 3.76 x 10" 11 (eV) 2 , 

bp 2 ~ 5.34 x 10~ 2 (eV) 2 . (57) 

So the physical scale k for 87 Rb becomes about 1.84 x 
10 2 nra. 



V. VORTEX SOLUTIONS IN 
GROSS-PITAEVSKII THEORY 



The two-component Gross-Pitaevskii theory is known 
to have non- Abelian vortices [13, HH • To obtain the vor- 
tex solutions in two-component Gross-Pitaevskii theory 
we first consider a straight vortex with the ansatz 



p = p(q), C 



fie) , ■ n- 

cos — — exp(— imp) 



sm • 



(58) 



5 10 15 20 

e 

FIG. 4: The untwisted vortex in the Gross-Pitaevskii theory 
of two-component BEC. Here we have put n = 1, and g is in 
the unit of k. Dashed and solid lines correspond to = 0.1 
and 0.2 respectively. 
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FIG. 5: The helical vortex in the Gross-Pitaevskii theory of 
two-component BEC. Here we have put m — n — 1, k = 
0.25/ k, and g is in the unit of k. Dashed and solid lines 
correspond to Sfi/fi = 0, and 0.1 respectively. 



With the ansatz is reduced to 



1 (\ - 9 n 2 ,n 2 „ 9N 9 f " 
P + -p-(-/ 2 + 7 -(^-V)sin^,, 

A 



1 



= 0. 



(59) 



Now, we choose the following ansatz for a helical vortex 

m 

/ fie) i ■ w 

I cos exp(— imp) \ 

P = Pie), C= /p) , (60) 

\ sm — — — exp(imkz) / 
and find that the equation l|38|l becomes 



9 ' -/' ( 7/ 2 + 2 



n 
? 



/ • + Q + 2^ / + - m 2 fc 2 - ^ 2 ) sin / 



0. 



(61) 



Notice that mathematically this equation becomes iden- 
tical to the equation of the straight vortex except 
that here 5/j 2 is replaced by 8[i 2 + m 2 k 2 . 



Now, with the boundary condition 

p(0) = 0, p(oo) = po, 

/(0) = 7T, /(oo) = 0, 



(62) 



we can solve l|61|l . With m — 0, n — 1 we obtain the 
straight (untwisted) vortex solution shown in Fig. but 



with m = n = 1 we obtain the twisted vortex solution 
shown in Fig. [5J Of course H38fl also admits the well- 
known Abelian vortices with £i = or (2 = 0. But ob- 
viously they are different from the non- Abelian vortices 
discussed here. 

The untwisted non- Abelian vortex solution has been 
discussed before 0, 0] , but the twisted vortex solution 
here is new |18| . Although they look very similar on the 
surface, they are quite different. First, when <5/i 2 = 0, 
there is no untwisted vortex solution because in this case 
the vortex size (the penetration length of the vorticity) 
becomes infinite. However, the helical vortex exists even 
when 5 /j, 2 = 0. This is because the twisting reduces the 
size of vortex tube. More importantly, they are physi- 
cally different. The untwisted vortex is made of a single 
vorticity flux, but the helical vortex is made of two vor- 
ticity fluxes linked together [Isj . 

In Skyrme theory the helical vortex is interpreted as 
a twisted magnetic vortex, whose flux is quantized due 
to the topological reason. The helical vortex in Gross- 
Pitaevskii theory is also topological, which can be viewed 
as a quantized vorticity flux |l8j. To see this notice that, 
up to the overall factor two, the potential C M introduced 
in i|39|l is nothing but the velocity potential (more 
precisely the momentum potential) of the doublet £ 0, 

m 

= -^(cos/ + 1)8^ - ^-(cos/ - l)d^z, (63) 
which generates the vorticity 



Huy = daV v - d v Vu 



f . 
= -sm 



fyiid^Qdutp - dugd^cp) 
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FIG. 6: The supercurrent i$ (in one period section in z- 
coordinate) and corresponding magnetic field Hg circulat- 
ing around the cylinder of radius g of the helical vortex 
in the Gross-Pitaevskii theory of two-component BEC. Here 
m = n = 1, k — 0.25/s, 5/i 2 — 0, and g is in the unit of k. 
The current density j$ is represented by the dotted line. 



FIG. 7: The supercurrent i$ and corresponding magnetic field 
H<p flowing through the disk of radius g of the helical vortex 
in the Gross-Pitaevskii theory of two-component BEC. Here 
m = n = 1, k — 0.25/k, 5fi 2 = 0, and g is in the unit of k. 
The current density jg is represented by the dotted line. 



+mk{d tl gd v z - d v gd^z)j . (64) 
This has two quantized vorticity fluxes, $ z along the z- 



$2 = / H^gdgdip = -2irn, 



(65) 



and <J>;p around the z-axis (in one period section from 
z = to z = 2n/k) 



- 1 - , 



Hzgdgdz = lum. 



(66) 



Clearly two fluxes are linked together. 



Furthermore, just as in Skyrme theory, these fluxes 
can be viewed to originate from the helical supercurrent 
which confines them with a built-in Meissner effect 



sin/ 



sin/' 



cos/ 



+mk(f+^Lf + -f)d,z 
smf g 



(67) 



This produces the supercurrents i$ (per one period sec- 
tion in ^-coordinate from z — to z — 2ir/k) around the 
z-axis 



2im sin / 
k g 



Q=0 



and i| along the z-axis 



iz = — 27rmfcgi/sin/ 



(68) 



(69) 



The vorticity fluxes and the corresponding supercurrents 
are shown in Fig. and Fig. [7| This shows that the 
helical vortex is made of two quantized vorticity fluxes, 
the $ z flux centered at the core and the flux sur- 
rounding it 01- This is almost identical to what we 
have in Skyrme theory. Indeed the remarkable similarity 
between Fig. and Fig. and Fig. 4 and Fig. 5 in 
Skyrme theory is unmistakable. This confirms that the 
helical vortex of two-component BEC is nothing but two 
quantized vorticity fluxes linked together. We emphasize 
that this interpretation holds even when the d/i 2 is not 
zero. 

The quantization of the vorticity l|65|l and (|66|l is 
due to the non-Abelian topology of the theory. To see 
this notice that the vorticity (|64(l is completely fixed by 
the non-linear sigma field h defined by £• Moreover, for 
the straight vortex n naturally defines a mapping ^(S* 2 ) 
from the compactified two-dimensional space S 2 to the 
target space S 2 . This means that the vortex in two- 
component BEC has exactly the same topological origin 
as the baby skyrmion in Skyrme theory. The only dif- 
ference is that the topological quantum number here can 
also be expressed by the doublet £ 



Qv = —^J e ij d i ^d j Cd 2 x = 



(70) 



Exactly the same topology assures the quantization of 
the twisted vorticity flux [l8j . This clarifies the topologi- 
cal origin of the non-Abelian vortices of Gross-Pitaevskii 
theory in two-component BEC. 

The helical vortex will become unstable unless the 
periodicity condition is enforced by hand. But just as in 
Skyrme theory we can make it a stable knot by making 
it a twisted vortex ring smoothly connecting two peri- 
odic ends. In this twisted vortex ring the periodicity 
condition of the helical vortex is automatically guaran- 



10 



teed, and the vortex ring becomes a stable knot. In this 
knot the n flux winds around m flux <I>^ of the he- 
lical vortex. Moreover the ansatz (|58|l tells that <f>£ is 
made of mainly the first component while is made of 
mainly the second component of two-component BEC. 
So physically the knot can be viewed as two vorticity 
fluxes linked together, the one made of the first compo- 
nent and the other made of the second component which 
surrounds it. 

As importantly the very twist which causes the insta- 
bility of the helical vortex now ensures the stability of 
the knot. This is so because dynamically the momentum 
mk along the z-axis created by the twist now generates 
a net angular momentum which provides the centrifugal 
repulsive force around the z-axis preventing the knot to 
collapse. 

Furthermore, this dynamical stability of the knot is 
now backed up by the topological stability. Again this 
is because the non-linear sigma field n, after forming a 
knot, defines a mapping ir 3 (S 2 ) from the compactified 
space S 3 to the target space S 2 . So the knot acquires 
a non-trivial topology ir 3 (S 2 ) whose quantum number is 
given by the Chern-Simon index of the velocity potential, 



Q = / ^kCd l C{d^d k C)d 3 x 



16tt 2 



( ,11 jiil'i- = run. 



(71) 



This is precisely the linking number of two vorticity 
fluxes, which is formally identical to the knot quantum 
number of Skyrme theory 0, 0, This assures the 

topological stability of the knot, because two fluxes linked 
together can not be disconnected by any smooth defor- 
mation of the field configuration. 

Similar knots in Gross-Pitaevskii theory of two- 
component BEC has been discussed in the literature 
[Til fl5l | . Our analysis here tells that the knot in Gross- 
Pitaevskii theory is a topological knot which can be 
viewed as a twisted vorticity flux ring linked together. 

As we have argued our knot should be stable, dy- 
namically as well as topologically. On the other hand 
the familiar scaling argument indicates that the knot in 
Gross-Pitaevskii theory of two-component BEC must be 
unstable. This has created a confusion on the stability of 
the knot in the literature 0G3- To clarify the confu- 
sion it is important to realize that the scaling argument 
breaks down when the system is constrained. In our case 
the Hamiltonian is constrained by the particle number 
conservation which allows us to circumvents the no-go 
theorem and have a stable knot 113, UM ■ 



VI. GAUGE THEORY OF TWO-COMPONENT 
BEC 



The above analysis tells that the non-Abelian vortex 
of the two-component Gross-Pitaevskii theory is noth- 
ing but a vorticity flux. And creating the vorticity 
costs energy. This implies that the Hamiltonian of 
two-component BEC must contain the contribution of 
the vorticity. This questions the wisdom of the Gross- 
Pitaevskii theory, because the Hamiltonian has no 
such interaction. To make up this shortcoming a gauge 
theory of two-component BEC which can naturally ac- 
commodate the vorticity interaction has been proposed 
recently In this section we discuss the gauge theory 
of two-component BEC in detail. 

Let us consider the following Lagrangian of U(l) 
gauge theory of two-component BEC [5j 



C = l 2 
X 

~2 



Di 



_7i_ 
2M 

l -H 2 

4 ^' 



(72) 



where = + igC^, and is the field strength of 
the potential C^. Two remarks are in order here. First, 
from now on we will assume 



Sn = 0, 



(73) 



since the symmetry breaking interaction can always be 
treated as a perturbation. With this the theory acquires 
a global U (2) symmetry as well as a local U (1) symmetry. 
Secondly, since we are primarily interested in the self- 
interacting (neutral) condensate, we treat the potential 
C M as a composite field of the condensate and identify 
with the velocity potential of the doublet £ 



c» = --c%c=-v li . 

9 9 



(74) 



With this the last term in the Lagrangian now repre- 
sents the vorticity l|64|l of the velocity potential that we 
discussed before 

Hp, = ~(d^d v C ~ d^d^C) = -H^. (75) 

This shows that the gauge theory of two-component 
BEC naturally accommodates the vorticity interaction, 
and the coupling constant g here describes the coupling 
strength of the vorticity interaction || . This vorticity in- 
teraction distinguishes the gauge theory from the Gross- 
Pitaevskii theory. 

At this point one might still wonder why one needs 
the vorticity in the Lagrangian l|72|) , because in ordinary 
(one-component) BEC one has no such interaction. The 
reason is that in ordinary BEC the vorticity is identi- 
cally zero, because there the velocity is given by the gra- 
dient of the phase of the complex condensate. Only a 
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non-Abelian (multi-component) BEC can have a non- 
vanishing vorticity. More importantly, it costs en ergy 
to create the vorticity in non-Abelian superfhrids jlflj . 
So it is natural that the two-component BEC (which is 
very similar to non-Abelian superfluids) has the vortic- 
ity interaction. Furthermore, here we can easily control 
the strength of the vorticity interaction with the coupling 
constant g. Indeed, if necessary, we could even turn off 
the vorticity interaction by putting g — oo. This justifies 
the presence of the vorticity interaction in the Hamilto- 
nian. 

Another important difference between this theory and 
the Gross-Pitaevskii theory is the U(l) gauge symmetry. 
Clearly the Lagrangian (|72|) retains the U(l) gauge in- 
variance, in spite of the fact that the gauge field is re- 
placed by the velocity field l|74|l . This has a deep impact. 
To see this notice that from the Lagrangian we have the 
following Hamiltonian in the static limit (again normal- 
izing to y/2M /%4>) 

n = \id,p) 2 + ^(m 2 - ic^ci 2 ) + \{ P 2 pD 2 



2// 

T' 



(76) 



Minimizing the Hamiltonian we have the following equa- 
tion of motion 

d 2 p-(\d i (\ 2 -\?d i e)p=±(p*-pi)p, 

{(d 2 - c f 5 2 C) + + -L;) A ;),C<>,; - d^diQ 

l V p gtp 2 - 

-C t ^c)(^-C t ^C)}c = o. (77) 



But factorizing £ by the U(l) phase 7 and CP 1 field £ 
with 



C = exp(i7)£, 



(78) 



we have 



\d,c\ 2 - ic^ci 2 = \d,e - i^i 2 

= -n ■ (dfj,n x dvh) = gH^, (79) 
so that we can rewrite 1)770 in terms of £ 

o 2 p- (m* - \^ds 2 )p = ^(p 2 - pDp, 

{(a 2 - £ta 2 o + + -Ljojo.cnx - d^m 

-tfditydi = 0. (80) 



Moreover we can express the Hamiltonian i|76[l com- 
pletely in terms of the non-linear sigma field n (or equiv- 
alently the CP 1 field £) and p as 

n= 1 -(d l p) 2 + 1 -p 2 (d l nf + ±(p 2 -plf 



^Pt{\{^ 2 + \p 2 {d l nf + \{p 



2 1)2 



+ 1^2^" X 



(81) 



This is because of the U(l) gauge symmetry. The U(l) 
gauge invariance of the Lagrangian l|72|l absorbs the U(l) 
phase 7 of (, so that the theory is completely described 
by £. In other words, the Abelian gauge invariance of 
effectively reduces the target space S 3 of £ to the gauge 
orbit space S 2 = S 3 /S 1 , which is identical to the target 
space of the CP 1 field £. And since mathematically £ is 
equivalent to the non-linear sigma field h, one can express 
(|76[) completely in terms of h. 

This tells that the equation l)77|) can also be expressed 
in terms of n. Indeed with H4Q(I . (|41|l . and (|79|l we can 
obtain the following equation from l|77|l [j| 

8^-1(8^= -(p'-pDp, 



p 



n x d 2 h + 2^^h x dih H — ^diHijdjfi 



gp 



= 0. 



(82) 



This, of course, is the equation of motion that one obtains 
by minimizing the Hamiltonian JSy. So we have two 
expressions, (17711 and l|82l) . which describe the equation 
of gauge theory of two-component BEC. 

The above analysis clearly shows that our theory of 
two-component BEC is closely related to the Skyrme the- 
ory. In fact, in the vacuum 



2 2 
P = Pa, 



1 



9 2 Po 



(83) 



the two equations l|llf) and l|82ll become identical. Fur- 
thermore, this tells that the equation Ij49(l of Gross- 
Pitaevskii theory is very similar to the above equation 
of gauge theory of two-component BEC. Indeed, when 
Sp 2 = 0, (15)l and ljg5jl become almost identical. This 
tells that, in spite of different dynamics, the two theories 
are very similar to each other. 



VII. TOPOLOGICAL OBJECTS IN GAUGE 
THEORY OF TWO-COMPONENT BEC 



Now we show that, just like the Skyrme theory, the 
theory admits monopole, vortex, and knot. We start 
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FIG. 8: The monopole solution in the gauge theory of two- 
component BEC. Here we have put A = 1 and r is in the unit 
of l/p 



FIG. 9: The non-Abelian vortex (dashed line) with m — 
0, n = 1 and the helical vortex (solid line) with m = n = 1 in 
the gauge theory of two-component BEC. Here we have put 
X/g 2 — 1, k — 0.64/k, and g is in the unit of k. 



from the monopole. Let 

-71 



p£ (7 = 0), 



(cos — exp(— iip) 
2 q 
sin - 



(84) With this the equation (|8Tj|l is reduced to 



and find 



n = £'a£ = r. 



(85) 



where (r, 0, <p) are the spherical coordinates. With this 
the second equation of (|80fl is automatically satisfied, and 
the first equation is reduced to 

2 1 A 

So with the boundary condition 

p(0) = 0, p(oo) = p , (87) 

we have a spherically symmetric solution shown in Fig. [S] 
Obviously this is a Wu-Yang type vorticity monopole 
dressed by the scalar field p ja, 0- In spite of the 
dressing, however, it has an infinite energy due to the 
singularity at the origin. 

Next we construct the vortex solutions. To do this we 
choose the ansatz in the cylindrical coordinates 

/ fie) i ■ w 

I cos exp(— imp) \ 

p = p(*), e=( . f f g) J, (88) 

\ sin — - — expyimkz) / 

from which we have 

(sin / cos (nip + mkz) \ 
sin / sin + mkz) , 
cos/ / 



(89) 



I. . 1( , 



P+-p-i[f 2 + (m 2 k 2 + —)sm 2 f)p 



t;(p 2 -Po>> 



g z P / 



■7 + 2- + (— + m fc ) — / 



1 „p ,T 

p p r 



TP 



In 2 

,2 



2 ; 2 

m k 



sin 2 / 

2 2 

.9 P 



(— + m k ) sin / cos / = 0. 



(90) 



Notice that the first equation is similar to what we have in 
Gross-Pitaevskii theory, but the second one is remarkably 
similar to the helical vortex equation in Skyrme theory. 
Now with the boundary condition 



P(fj) = 0, p(oo) = po, 
/(D) = tt, f(oo) = 0, 



(91) 



we obtain the non-Abelian vortex solution shown in 
Fig. El 

The solution is similar to the one we have in Gross- 
Pitaevskii theory. First, when m = 0, the solution de- 
scribes the straight non-Abelian vortex. But when m is 
not zero, it describes a helical vortex which is periodic 
in z-coordinate 5]. In this case, the vortex has a non- 
vanishing velocity current (not only around the vortex 
but also) along the z-axis. Secondly, the doublet £ starts 
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FIG. 10: The supercurrent i$ (in one period section in z- 
coordinate) and corresponding magnetic field Hg circulating 
around the cylinder of radius g of the helical vortex in the 
gauge theory of two-component BEC. Here m = n = 1, 
\/g 2 — 1, k — 0.64/k, and g is in the unit of k. The current 
density j$ is represented by the dotted line. 
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FIG. 11: The supercurrent i$ and corresponding magnetic 
field Htp flowing through the disk of radius g of the helical 
vortex in the gauge theory of two-component BEC. Here m = 
n — 1, X/g 2 = 1, k = 0.64/k, and g is in the unit of k. The 
current density j$ is represented by the dotted line. 



from the second component at the core, but the first com- 
ponent takes over completely at the infinity. This is due 
to the boundary condition /(0) = ir and /(oo) = 0, which 
assures that our solution describes a genuine non-Abelian 
vortex. This confirms that the vortex solution is almost 
identical to what we have in the Gross-Pitaevskii theory 
shown in Fig. 5. All the qualitative features are exactly 
the same. This implies that physically the gauge the- 
ory of two-component BEC is very similar to the Gross- 
Pitaevskii theory, in spite of the obvious dynamical dif- 
ferences. 

Clearly our vortex has the same topological origin as 
the vortex in Gross-Pitaevskii theory. This tells that, 
just as in Gross-Pitaevskii theory, the non-Abelian heli- 
cal vortex here is nothing but the twisted vorticity flux of 
the CP 1 field £ confined along the z-axis by the velocity 
current, whose flux is quantized due to the topological 
reason. The only difference here is the profile of the vor- 
ticity, which is slightly different from that of the Gross- 
Pitaevskii theory. Indeed the solution has the following 
vorticity 

1 ^ 



H, = -H, 



Q 2gg 



H{p — —Hg Z 



mk 
^9 



/sin/, 



(92) 



which gives two quantized vorticity fluxes, a flux along 
the z-axis 

4>z = J H ev dgdip 

= -— [ (0 c ft a e _ dip t%t)dQ = -—, (93) 

9 J 9 
and a flux around the z-axis (in one period section from 



to 2n/k in z-coordinate) 



2m 
9 



Hgz dgdz 



(d e ?d z z-d z tfd e t)=i = 



dg 2itm 



9 



(94) 



This tells that the vorticity fluxes are quantized in the 
unit of 2ir/g. 

Just like the Gross-Pitaevskii theory the theory has a 
built-in Meissner effect which confines the vorticity flux. 
The current which confines the flux is given by 



sin/ 



'Ml 



3 fi ^yH^y 

sin/ g 



ik(f+ C -^.f 2 + -f)d,z 
sin/ Q 

d^jn = 0. 



(95) 



This produces the supercurrents (per one period sec- 
tion in z-coordinate from z = to z = 2n/k) around the 
z-axis 



7rn sin / 
gk g 



g=0 



and i£ along the z-axis 



mk ■ . 

-t QJ sin/ 

9 



(96) 



(97) 



Q=0 



The helical vorticity fields and supercurrents are shown 
in Fig.llOland Fig. ^2 The remarkable similarity between 
these and those in Skyrme theory (Fig. [21 and Fig.[2|) and 
Gross-Pitaevskii theory (Fig. Eland Fig. 0) is unmistak- 
able. 
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With the ansatz l|89[l the energy (per one periodic 
section) of the helical vortex is given by 

4tt 2 f°° fl „ 1 „/, 



1 ,n 



g2p2 g2 



+m 2 k 2 ) sin 2 . f)f 2 + ('- 



+m 2 k 2 )sin 2 f)+^(p 2 -p 2 ) 2 }gdg 



J2 poo 



If dp, 2 1 



+4^(4+™ 2 ^w/)(fr 

5/3 ar ax 



+ m 2 K 2 k 2 ) sin 2 / 



+ -( / 5 2 - l) 2 ^dx. 



(98) 



One could calculate the energy of the helical vortex nu- 
merically. With m = n = 1 and fc = 0.64 k we find that 
the energy in one period section of the helical vortex in 
87 Rb is given by 

E ~ 51 ~ 4.5 x 1CT 10 ey 

~ 0.7MHz, (99) 
which will have an important meaning later. 



VIII. VORTICITY KNOT IN 
TWO-COMPONENT BEC 



The existence of the helical vortex predicts the exis- 
tence of a topological knot in the gauge theory of two- 
component BEC, for exactly the same reason that the 
helical vortices in Skyrme theory and Gross-Pitaevskii 
theory assure the existence of knots in these theories. To 
demonstrate the existence of knot in the gauge theory 
of two-component BEC we introduce the toroidal coor- 
dinates (77,7, p) defined by 

x = — sinh i] cos <p, y — — sinh 77 sin ip, 
a 

z = — sin 7, 



D = cosh rj — cos 7, 



ds 2 = -^{drj 2 + d-y 2 + sinh 2 rjdip'' 



d 3 x — — j sinh ndrid'y dip, 
D A 



(100) 



where a is the radius of the knot defined by 77 = 00. 
Notice that in toroidal coordinates, 77 = 7 = represents 
spatial infinity of R 3 , and r\ = 00 describes the torus 
center. 



Now we choose the following ansatz 



V2 



exp(— inui(r), 7)) 



sm — - — e~xjp(imip) 



(101) 



With this we have the velocity potential 
Cn = - — (cos / - l)9 M y> 



2ff 



2.9 

(cos / + l)<9 M w, 



(102) 



which generates the vorticity 

Hfj, v = d^C v ~ d v C^, 
H V1 — —Ksinf, H lv = — sin fd 1 f, 



2.9 



25 



H<pri = - ■z- sin fd v f 
^9 



(103) 



where 



K = d^fd^oj - d 1 fd v uj. (104) 
Notice that, in the orthonormal frame (17,7, ip), we have 

Cfi = -77— ( cos / + i)dv u > 

2ga 
nD 



Cj = — - — (cos / + l)d y uj, 



6 



2ga 

mD 
2ga sinh 77 



(cos / - 1), 



(105) 



and 



77 D 2 

II,, ^7 A MU/. 



711 ~ 2ga 2 

lv 25a 2 sinh 77 
- mD 2 . 

^ = -2^sinhr ? Sln/9 " / - 



sin fdjf, 



(106) 



From the ansatz (llOlfl we have the following equations 
of motion 



* sinh 77 



D 



1 r 



(dr,f) 2 + (ft,/) 5 



+ (n 2 ((d, 1 uj) 2 + (d J u J ) 2 ) + - 
v si 



m 



sinh 2 77 



D 



sin 2 / 



Aa 2 / ., 



Q 2 + d 2 + I - 
" 7 V sinh 77 £> 

-(n 2 ((^) 2 + (^) 2 ) + 



2L> 2 

^ cosh 77 sinh 77 ^ sin 7 



9. 



2 I sin / cos / 
sinh rj ' 
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FIG. 12: (Color online). The p profile of BEC knot with m FIG. 13: (Color online). The / profile of BEC knot with 
= n = 1. Here we have put X/g 2 = 1, and the scale of the m — n = 1. Here we have put X/g 2 = 1. 
radius a is k. 
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(107) 



where 



A = n 2 K 2 



sinh 77 



(W) 2 + (a 7 /) ; 



+n 2 K 



^, cosh 77 ^ sinh 77.^ 
sinh 77 D 



sin 7 
—p-d v uj 



sinh 77 



cosh 77 sinh 77,^ 
sinh 77 D 



d 2 + d 2 
sin 7 



a. 



/, 



C = ^if9 7 / - dr,fd 7 K 
sin 7 



cosh 77 sinh 77 
( ^inh^ + _ D~ )d7 



Since 77 = 7 = represents spatial infinity of R 3 and 
77 = 00 describes the torus center, we can impose the 
following boundary condition 



p(0,0)=p , 0(00,7) =0, 
/(0 )7 ) = 0, /(oo )7 )=7r, 
w(»7,0)=0, a;(7 ? ,27r) = 2 7 r ) 



(108) 



to obtain the desired knot. From the ansatz (|101|l and 
the boundary condition l|108(l we can calculate the knot 
quantum number 



Q = ^2 / K sin fdr)djdtp 



An 



sin fdfdiu 



(109) 



where the last equality comes from the boundary condi- 
tion. This tells that our ansatz describes the correct knot 
topology. 

Of course, an exact solution of (|107(l with the bound- 
ary conditions (|108J) is extremely difficult to find 0, E| • 
But here we can obtain the knot profile of p, /, and ui 
which minimizes the energy numerically. We find that, 
for m = 71 = 1, the radius of knot which minimizes the 
energy is given by 



1.6k. 



(110) 



From this we obtain the following solution of the light- 
est axially symmetric knot in the gauge theory of two- 
component BEC (with m = n = 1) shown in Fig. 1121 
Fig. and Fig. [21 With this we can obtain a three- 
dimensional energy profile of the lightest knot (Unfor- 
tunately we can not show the profile here because the 
volume of the eps-file is too big) . 

We can calculate the vorticity flux of the knot. Since 
the flux is helical, we have two fluxes, the flux $ 7 passing 
through the knot disk of radius a in the xy-plane and the 
flux which surrounds it. From (|106J) we have 



$ 7 = 



- a 2 sinh 77 



H. 



D 2 



= -— sin fd v fdi]dip = - 

^9 J~i = -K 



drjdkp 
27rm 



a 



(111) 
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FIG. 14: (Color online). The uj profile of BEC knot with 
m = n = 1. Notice that here we have plotted to — 7. Here 
again we have put X/g 2 — 1. 
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FIG. 15: (Color online). The Q-dependence of axially sym- 
metric knot. The solid line corresponds to the function 
EoQ 3 ^ with Eo = E(l, 1), and the red dots represent the 
energy E(m, n) with the different Q — mn. 



and 



(112) 



n I TS • (A A 2lTn 

= — K sin / and'y = . 

2g J g 

This confirms that the flux is quantized in the unit of 
2n/g. As importantly this tells that the two fluxes are 
linked, whose linking number is fixed by the knot quan- 
tum number. 

Just as in Gross-Pitaevskii theory the vorticity flux 
here is generated by the helical vorticity current which is 
conserved 



nD 2 



3n = 
nD 2 

mD 2 
~2g~a^ 



sm 



2ga 2 
f(dn 



in /(8 



sin 7 



^ D 
cosh r\ sinh i] 



sinh 77 D 
cosh r\ sinh ?y 
sinh 7/ D 



sin fdnf 



D J - 
drfp = 0. 



(113) 



Clearly this supercurrent generates a Meissner effect 
which confines the vorticity flux. 

From 176( 1 and (|1(J1(I we have the following Hamilto- 
nian for the knot 



(9 7 p) 2 



(9,/) 2 + (9 7 /) 2 + (n 2 ((a,c) 2 + (c^) 2 ) 



sinh 77 



sm 

-D 4 
~Sg 2 a' 



-A sin /. 



(114) 



With this the energy of the knot is given by 

„3 



E = / TL — - smhridndjclip 
J D 6 

- J ^^753 sinhridrjdjdip, 



(115) 



where 



n 



K 2 D< 

2a 2 



{(d n p) 2 + (d,pf 



(d n f) 2 + (d 7 f) 2 + (n 2 ((d n co) 2 + (d^) 2 ) 
in 2 /} + i(/5 2 -l) 2 



sinh 2 r\ 



A k 4 D 4 



Asm 2 f. 



(116) 



Minimizing the energy we reproduce the knot equation 
(tTU7l) . 

From this we can estimate the energy of the axially 
symmetric knots. For the lightest knot (with m = n = 1) 
we find the following energy 



£-54 



Po 



Va 

~ 0.75MHz 



4.8 x l(T 10 ey 



(117) 



One should compare this energy with the energy of the 
helical vortex l|99(l . Notice that the lightest knot has 
the radius r ~ 1.6k. In our picture this knot can be 
constructed bending a helical vortex with k ~ 0.64/k. 
So we expect that the energy of the lightest knot should 
be comparable to the energy of the helical vortex with 
k ~ 0.64/k. And we have already estimated the energy 
of the helical vortex with k ~ 0.64/k in The fact 

that two energies are of the same order assures that the 
knot can indeed be viewed as a twisted vorticity flux ring. 
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As we have remarked the Q-dependence of the energy 
of Faddeev-Niemi knot is proportional to Q 3 / 4 0, Il3j| . 
An interesting question is whether we can have a simi- 
lar Q-dependence of energy for the knots in BEC. With 
our ansatz we have estimated the energy of knot numer- 
ically for different m and n up to Q = 6. The result is 
summarized in Fig. 1151 which clearly tells that the en- 
ergy depends crucially on m and n. Our result suggests 
that, for the minimum energy knots, we have a similar 
(sub- linear) Q-dependence of energy for the knots in two- 
component BEC. It would be very interesting to establish 
such Q-dependence of energy mathematically. 



IX. DISCUSSION 

In this paper we have discussed two competing theo- 
ries of two-component BEC, the popular Gross-Pitaevskii 
theory and the U(l) gauge theory which has the vor- 
ticity interaction. Although dynamically different two 
theories have remarkably similar topological objects, the 
helical vortex and the knots, which have a non-trivial 
non-Abelian topology. 

We have shown that the U(l) x U(l) symmetry of 
two-component BEC can be viewed as a broken U(2) 
symmetry. This allows us to interpret the vortex and 
knot in two-component BEC as non-Abelian topological 
objects. Furthermore, we have shown that these topolog- 
ical objects are the vorticity vortex and vorticity knot. 

A major difference between the Gross-Pitaevskii the- 
ory and the gauge theory is the vorticity interaction. In 
spite of the fact that the vorticity plays an important role 
in two-component BEC, the Gross-Pitaevskii theory has 
no vorticity interaction. In comparison, the gauge the- 
ory of two-component BEC naturally accommodates the 
vorticity interaction in the Lagrangian. This makes the 
theory very similar to Skyrme theory. More significantly, 
the explicit U{1) gauge symmetry makes it very simi- 
lar to the theory of two-gap superconductor. The only 
difference is that the two-component BEC is a neutral 
system which is not charged, so that the gauge interac- 
tion has to be an induced interaction. On the other hand 
the two-gap superconductor is made of charged conden- 
sates, so that it has a real (independent) electromagnetic 
interaction [20| . 

As importantly the gauge theory of two-component 
BEC, with the vorticity interaction, could play an im- 
portant role in describing multi-component superfluids 

UM- In fact we believe that the theory could well 
describe both non-Abelian BEC and non-Abelian super- 
fluids. 

In this paper we have constructed a numerical solu- 
tion of knot in the gauge theory of two-component BEC. 
Our result confirms that it can be identified as a vortex 
ring made of a helical vorticity vortex. Moreover our re- 



sult tells that the knot can be viewed as two quantized 
vorticity fluxes linked together, whose linking number be- 
comes the knot quantum number. This makes the knot 
very similar to Faddeev-Niemi knot in Skyrme theory. 

We close with the following remarks: 

1. Recently a number of authors have also established the 
existence of knot identified as the "skyrmions" in Gross- 
Pitaevskii theory of two-component BEC [Til Il5| , which 
we believe is identical to our knot in Gross-Pitaevskii the- 
ory. In this paper we have clarified the physical meaning 
of the knot. The knot in Gross-Pitaevskii theory is also 
of topological origin. Moreover, it can be identified as a 
vorticity knot, a twisted vorticity flux ring, in spite of the 
fact that the Gross-Pitaevskii Lagrangian has neither the 
velocity C M nor the vorticity _ff M „ which can be related 
to the knot. 

2. Our analysis tells that at the center of the topologi- 
cal vortex and knot in two-component BEC lies the baby 
skyrmion and the Faddeev-Niemi knot. In fact they are 
the prototype of the non-Abelian topological objects that 
one can repeatedly encounter in physics [H, El EH- 
This suggests that the Skyrme theory could also play 
an important role in condensed matter physics. Ever 
since Skyrme proposed his theory, the Skyrme theory 
has always been associated to nuclear and/or high en- 
ergy physics. This has lead people to believe that the 
topological objects in Skyrme theory can only be real- 
ized at high energy, at the GeV scale. But our analysis 
opens up a new possibility for us to construct them in 
a completely different environment, at the eV scale, in 
two-component BEC 0.0- This is really remarkable. 

3. From our analysis there should be no doubt that 
the non-Abelian vortices and knots must exist in two- 
component BEC. If so, the challenge now is to ver- 
ify the existence of these topological objects experimen- 
tally. Constructing the knots might not be a simple 
task at present moment. But the construction of the 
non-Abelian vortices could be rather straightforward, 
and might have already been done [j, |22| • Identifying 
them, however, may be a tricky business because the 
two-component BEC can admit both the Abelian and 
non-Abelian vortices. To identify them, one has to keep 
the following in mind. First, the non-Abelian vortices 
must have a non-trivial profile of f(g). This is a crucial 
point which distinguishes them from the Abelian vor- 
tices. Secondly, the energy of the non-Abelian vortices 
must be bigger than that of the Abelian counterparts, 
again because they have extra energy coming from the 
non-trivial profile of /. With this in mind, one should 
be able to construct the non-Abelian vortices in the new 
condensates without much difficulty. We strongly urge 
the experimentalists to meet the challenge. 
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